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ABSTRACT 

We investigate the feasibility of extracting the gravitational nanolensing signal due to the presence 
of subsolar mass halos within galaxy-sized dark matter halos. We show that subsolar mass halos in 
a lensing galaxy can cause strong nanolensing events with shorter durations and smaller amplitudes 
than microlensing events caused by stars. We develop techniques that can be used in future surveys 
such as Pan-STARRS, LSST and OMEGA to search for the nanolensing signal from subsolar mass 
halos. 

Subject headings: galaxies: halos - gravitational lensing - theory: dark matter 



1. INTRODUCTION 

In cold dark matter (CDM) cosmologies, dark matter 
halos form in a hierarchical manner. Small dark mat- 
ter halos form first, and their subsequent merger results 
in the formation of larger systems. Dark matter halos 
contain a remnant of this process in the form of sub- 
structure, the self-bound smaller systems that survive 
the hierarchical assembly of their host halo. 

The scale of the first objects that form in the Universe 
is set by the nature of the dark matter particle. An 
observation of a cutoff scale in the dark matter power 
spectrum (either in the present-day Universe or in ear- 
lier epochs) would provide an insight into the particle 
physics properties of the dark matter, as well as infor- 
mation about the dynamics and evolution of non-linear 
systems at the dawn of structure formation. 

On scales greater than ~ lO^M©, luminous galaxies 
form in the center of the dark matter potential well 
(e.g., dwarf galaxies, galaxies, and clusters of galaxies). 
On scales less than ^ 10^ Mq the collisional nature of 
baryons prevents the formation of visible objects, making 
the use of luminous matte r as a tracer of the dark mat- 
ter distribution ineffective (Strigari et al. 12008 Kravtsov 



2010| [Strigari et al.||2010| ). I'hus, the existence ot selt- 
bound dark matter systems with masses less than ap- 
proximately 10^ Mq is very difficult to establish. This 
presents a very serious problem in attempts to under- 
stand the small-scale structure of dark matter halos. 

Weakly Interacting Massive Particles (WIMPs) are 
well established theoretically and offer a natural dark 



matter re lic density contribution (Jungma n et al.| 1996 
Bertone et al. 2005) . They arise in theories which in- 



troduce physics beyond the Standard Model and have 
masses of order of tens of GeV and an interaction 
strength with Standard Model particles of the order of 
the weak interaction. The free streaming of WIMPs af- 
ter kinetic decoupling in the early Universe (momentum- 
changing interactions) introduces a cutoff in the matter 
power spectrum typically at a sub-solar mass scale be 
low ^ IQ-'^Mq dSchmid et al.||l999[ [Green et al.||2004 



2005). This cutott scale has been theoretically explorec 
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in the particular case of supersymmetric WIMPs (Hof- 



mann et al.|[ 2QQT^ 'Chen et al.|[2QQT| |Profumo et al.. ^ 

Loeb Za ldarriaga 2005; Martinez et al.||2009| ). While 
initially the first dark matter halos had masses near the 
cutoff scale, it is unclear whether these objects survived 
the hierarchical growth of their host halo ([Die mand et al. 
^QQ5l|Berezinsky et al.| 2 0031 [Green & Goodwin 2007, 



erdt et al.|2007 Ishiyama et al.|2010p . Numerical siniula- 
tions, as well as sophisticated analytical arguments, can- 
not at present address this issue conclusively, and there- 
fore any potential observational signature of the presence 
of subsolar mass halos is valuable. 

The possibility of the detection of subsolar mass dark 
matter halos by indirect and direct detection experiments 
has been studied extensively in literature. Indirect de- 
tection experiments search for the products of dark mat- 
ter annihilation. As this process is proportional to the 
integral of the square of the dark matter density dis- 
tribution over volume, high density regions enhance the 
annihilation rate relative to low density regions. Sub- 
solar mass dark matter halos have very high densities 
as they are formed at extremely high redshifts (their 
densities reflect the mean dark matter density of the 
Universe at the redshift of formation). The probabil- 
ity of detecting the indirect detection effects of subsolar 
mass subhalos present in present day dark matter ha- 
los has been studied in the context of individual detec- 



tion ( Diemand et a l.|2005' Pieri et al.|2005 [Koushiap pas[ 
20061) , as well as their contribution to the diffuse ga mma- 



ray background ([Ando et al.[[2008 
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2010) 



Kamionk owski et al. 
2009; Baxter et al. 



PierTet af.[[:^00^[ [Lee et a. . 

In addition, a subject that received considerable 
attention is the role of subsolar mass substructure in set- 
ting the amplitude of the annihilation flux i n halos, the 
so-called "boost" factor (Bergstrom et al. 1999, 1998 
Calcaneo-Roldan fc Moore, 2000.; Be rezinsky et al. 2(17^ 
Stoehr et al.|2003|[Baltz et al.|20'07[ [D ieman d et al.|2007 



Berezinsky et al.''2007^ » Strigari'eFaT.||2007i |Kuhlen et al. 
2008, 2009; Kamionkow ski et al.[[201Q| ^ 

On the other hand, direct detection experiments are 
minimally sensitive to the presence of subsolar mass ha- 
los in the Milky Way. Inspired by a nalytical arguments 
(Kamionkowski & Koushiappas 2008^, numerical simu- 
lations show that the direct detection experiments are 
not sensitive to subsolar mass subhalos for two reasons 
( [Vogelsberger et al.,^2009) . First, the volume occupied 
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by subsolar mass substructure at the present epoch is 
extremely smah (thus the probabihty of an interaction 
with a subsolar mass halo during the course of a direct 
detection experiment is negligible). Second, the smooth 
dark matter distribution does not retain any informa- 
tion about the presence of subsolar mass halos (even as 
tidal streams) due to the very efficient mixing of tidally 
stripped material since the dynamical timescale at the 
solar radius i s a very small fraction of the age of the 
Milky Way ( Vogelsberger fc White||2QlQD . 

Finally, energetic neutrinos from the Sun and the Earth 
could potentially hold a signature of their interactions 
with subsolar mass halos along the Sun's Galactic or- 
bit: the signal s amples the past history of the local d ark 



matter density ( Koushiappas fc Kamionkowski|[2QQ9 ) . A 
past interaction with a subsolar mass halo could give rise 
to an elevated signal of energetic neutrinos, though such 
an effect is sensitive to the equilibration timescale of the 
dark matter particle, the survival of subsolar m ass halos 
in the Milky Way, and their internal structure (iKoushi-l 
appas fc Kamionkowski||2QQ9[ Serpico fc Bertone 2010). 



Gravitational lensing otters a new possibility for the de 
tection of subsolar mass dark matter halos. Light from 
distant objects is deflected by matter along the line of 
sight. In some cases, such as a sufficiently large inter- 
vening galaxy, the observer will see multiple images of 
the same source; this is referred to as 'strong' gravita- 
tional lensing. Strong lens systems with quasar sources 
and galaxy-sized lenses have be en shown to be sensi- 
tive to stars (Chang fc Refsdalll 19791 |Paczynski||1986 



tive to stars (U nang &^ neisdal [ibiYbi ■raczynsKi||lbioc 
Kayser et al.||1986| [Schneider <fe Weissl|1987[ |lrwin et a. 



Paczyn- 



,19891 IWambs^anss et al.||1990| [Wambsganss 

ski||199ip , as well as satellite g alaxy-sized dark matter 
subhalos in the lensing galaxy jMao fc Schneider" 1998 
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Metcalf fc Madau 2001; Metcall <fe Zhao||2002 ; Dalai 
Koch anek 2002; Chiba 2002; Kochan ^QrHal al 2004). 
I'he sensitivity of lensing systems to the larger dark mat- 
ter subhalos is a well-known prob e of the nature of the 
dark matter particle (see, e.g., |Moustakas et al. 2009j). 
These systems could also be sensitive to much sma ller 
dark matter substructure (Lewis fc Gil-Merino 2006). If 
so, strong gravitational lensing might provide the best 
estimate for detecting subsolar mass DM halos since the 
lensing signal is sensitive only to mass and does not de- 
pend upon the nature of the dark matter particle. 

In this paper, we focus on this new possibility of us- 
ing gravitational lensing as a probe of the existence of 
subsolar mass subhalos in dark matter halos. In Sec. [2] 
we present an order of magnitude estimate of the effect 
based on a simple toy model. As the main source of con- 
fusion in such an investigation is the presence of stars 
in the lens, we construct a mock lens system and inves- 
tigate the lensing effects of stars and dark matter halos 
in Sec. [3l We then create mock observations in Sec. [4] 
and search for the signatures of subsolar mass dark mat- 
ter halos in Sec. [5] In Sec. [6] we discuss challenges to 
measuring subsolar mass dark matter halos in real ob- 
servations, and we discuss the prospects for detection in 
current and future data sets in Sec. 

2. SIMPLE ESTIMATES 

In strong gravitational lens systems, a distant source 
is multiply-imaged by an intervening object. In cases for 
which the source is a quasar and the lens is a galaxy. 



small sub-galactic objects can cause large perturbations 
to the magnification of individual images. Investigating 
such effects is computationally intensive, so here we dis- 
cuss in a simplified way the prospects for finding such 
large magnification perturbations due to subsolar mass 
dark matter halos. 

A simple proxy for large magnification perturbations 
is the splitting of an image into subimages. The total 
magnification of the subimages is at least as big as the 
magnification in the absence of image splitting. In addi- 
tion, multiple subimages imply the existence of solutions 
of the lens equation for which the magnification of re- 
sulting subimages can be arbitrarily high - i.e., critical 
points. 

When the perturbers in the lens are stars, multi- 
ple subimages are always present if the star and the 
source are sufficiently aligned because stars are point 
masses. Stars cause image splittings at roughly the 
micro-arcsecond (/ias) scale, giving the effect the name 
'microlensing.J^ Subsolar mass dark matter halos, how- 
ever, are not point masses. A sufficient condition for 
multiple subimages is that the surface mass density is 
greater than the critical density of lensing. 



^total — S/^crit ^ 1- 



(1) 



Here we specify /^totai to emphasize the fact that can 
have contributions from many sources, such as subsolar 
mass halos or a smooth dark matter component. 

Let us consider some toy parameters for a subsolar 
mass halo. At the redshift of collapse, a ~1 so lar mass 
{Mq) halohas a virial radius of ~1 parsec (pc) ( [Koushi- 
appas| [20091 ) • for mation, studies suggest that the 
halo is described by a Navarro et al. (1997, NFW) pro- 
file, p(r) oc l/[x(l + x)^] where x = r/rg and r is the 
three-dimensional radius an d r,s is the scale radius ( |Die-| 
mand et al]|2005[ [20061 ). [Ishiyama et al.[ ( |2010[ ) tmd 
p oc r~^"^ using higher resolution simulations, but ignor- 
ing the larger scale fluctuations that could allow halos 
to form at earlier redshift and with higher central den- 
sities. Let us approximate the projected mass profile of 
the halo as a power law, S(i?) oc where R is the 

projected radius and a > 0. We consider a lensing galaxy 
at z = 0.3 and source quasar at 2: = 2. 

At the position of the images in strong lens systems, 
^totai ~ 0.5. If all the mass is in solar mass halos, then 
^ 3000 halos will overlap at any given point. Now let us 
separate the contribution of one halo from all the other 
halos so that A^totai = ^ihaio + ^back- Givcu the high av- 
erage number density of halos, the average contribution 
of all but one halo at a point is hZhack ~ 0.5. If we then 
consider the surface density contribution of a single halo, 
there is an area very close to the center of this halo where 
the surface density is /^ihaio > 0.5. This occurs at a ra- 
dius of 0.002, 0.0002, and 7 x 10"^ pc for a halo profile 
slope of a = 1.5, 1, and 0.5 respectively. An estimate 
of the probability of a source point falling into the area 
of /^totai > 1 is 0.03, 3 X 10-"^, and 4 x 10"^^ for a halo 
profile slope of a = 1.5, 1, and 0.5 respectively. 

The area of /^ihaio ^ 0.5 also corresponds roughly to 
the size of the deflection caused by the subsolar mass 
dark matter halos and is ~ 0.1 /ias for a = 1.5. The mass 



For a comprehensive discussion of microlensing, see |Wambs^ 
|ganss| ( [2006| >. 
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function of subsolar masses extends significantly below a 
solar mass, and the typical deflections for smaller halos 
will be correspondingly smaller and encompass the nano- 
arcsecond (nas) scale. As the lensing effects of stars are 
generically referred to as 'microlensing,' we refer to the 
effects of subsolar mass dark matter halos as 'nanolens- 
ing.' 

In this estimate, we made a significant leap by assum- 
ing /^back = 0.5, when in reality it has significant vari- 
ance. For a conservative estimate, we can set the back- 
ground to A^back = and consider the area around a sin- 
gle halo where A^ihaio ^ 1- Then the probabilities drop to 
0.01, 7 X 10-^ and 3 x 10-^^ for slopes of a = 1.5, 1, and 
0.5 respectively. Clearly, for a steep halo profile (a > 1) 
and a sufficient amount of observations, the likelihood 
of observing a significant nanolensing event is not small. 
Subsolar mass dark matter halos with power-law slopes 
shallower than a ~ 1, on the other hand, are unlikely to 
produce any observable effects. 

3. MICRO- AND NANOLENSING SIMULATIONS 

We create a mock lens system similar to typical ob- 
served systems. The lens galaxy is at z = 0.3 and the 
quasar source is at redshift z = 2. The magnification of 
each image in a smooth lens system is given by 



Mock Image 


hi 


7 


/i 


Minimum (M) 


0.475 


0.425 


10.5 


Saddle (S) 


0.525 


0.575 


-9.5 



(1 



■r 



(2) 



where is the local projected mass density and 7 is the 
external shear. Images in which l — Hc — ^y > are minima 
in the time travel surface and are positive parity images, 
while those with 1 — tz — j < are saddle points and are 
negative parity images. We create a minimum (M) and a 
saddle point (S) image that in the case for which all the 
mass is distributed in a smooth component (no stars, no 
subsolar mass DM halos) have a magnification typical for 
lens system with four images, /i ^ 10. The parameters 
for the mock images are chose n to be the same as in 
iSchechter & Wambsganssl (|2002|) and are shown in Table 

^ — 

Images in a four image system are located at ~3% 
of the projected virial radius of the lensing halo, a dis- 
tance where a significant fraction of the projected mass 
is expected to be in stars. Studies have found a range 
of values for k^Jk. For example, the study of Koop- 



mans et al. (2006) of the Sloan Lens ACS Survey found 
an average stellar mass fraction of 0.75 inside of the 
Einstein radius. Studies of lens systems with quasar 
sources have found gen erally lower stella r mass fractions: 
> 0.5 (Q2237+0 305; iKochanek l [2QQ4| , 0.08-0.15 (PG 
1115+080; Morgan et al.''200 8|/ (11^X3 (PG 1115+08 
and SDSS 0924+0219; Schech teTfc Wamb sgajiss| |2004D , 
- 0.1 (PG 1115+080; Pooley et al. 2009), 0.05 (|Medi- 



avilla et al. [2009 ) for a sample o f lens systems, and U.l 
(RXJ 1131-l!j31; ]Dai et al.||2QTQ ). 



We put half of the projected mass in each mock image 
realization in stars, hc^ = 0.5/^, a relatively high value and 
potentially a large contaminant to the signal from subso- 
lar mass dark matter halos. Stars are modeled as point 
mass perturbers with a Salpeter mass function ( |Salpeter| 



1955), dn/dm ~ m~^-'^^ with mass limits of O.OIM© to 
IMq. Stars are placed at random within an area with 
radius ~ 20 times greater than the Einstein radius of a 
solar mass star. 



TABLE 1 
Simulation parameters 



We introduce subsolar mass halos in the simulations in 
the following way. If subsolar mass dark matter halos are 
resistant to tidal disruption, they could make up the ma- 
jority of the dark matter halo of the Galaxy. The effects 
of tidal disruption on subsolar mass halos are still under 
debate, though recent studies suggest that the central 
cores will survive all but the densest regions (< 20pc 
from the G alactic center) of the Milky Way ( Ishiyama] 



et al. 2010). Based on this result and for the sake ol 
simplicity, we put the remaining half of the projected 
mass into subsolar mass halos. We model subsolar mass 
dark matter halos with power-law density profiles where 
^{R) ^ R~^-^. As discussed earlier, this may be steeper 
than expected and makes them more effective perturbers. 
A solar mass dark matter halo has a radius of 1 par sec. 
For smaller halos, the radius is scaled with the mass, 
M (X r^. The subsolar mass halo mass function is mod- 
eled as dn/dm ~ with mass range of lO^^M© to 
IMq. These halos are distributed randomly within the 
same area as stars. While stars, as point masses, al- 
ways induce microlensing as long as light from the source 
passes close enough to the star, subsolar mass halos may 
not act alone. Each realization has significantly more 
subhalos than stars (~ 10^ compared to ^ 10^), and 
many halos overlap. As a result, the mass density from 
every halo adds up and microlensing events are induced 
where the combined mass density of many subhalos ex- 
ceeds the critical density for lensing. 

We generate seven realizations for each mock saddle 
image (labeled Shi to Sh7), and five realizations for each 
mock minimum image (Mhl to Mh5). For comparison, 
we also generate the same number of realizations contain- 
ing only stars but no subsolar mass dark matter halos - 
in this case all the dark matter is in a smooth component 
(SI to S7 and Ml to M5). 

We perform a lensing simulation using an inverse ray- 
tracing code. The position of the source, /3, and the 
images, 6>, are related by the lens equation, j3 = — 
V(j){6)^ where (j) is the lensing potential. In the case of 
circular symmetry in the lens, V(j) oc M, where M is the 
mass of the lens. As the lens equation is mult i- valued 
and several image positions may correspond to a single 
source position, it is easiest to solve by establishing a grid 
of image positions and calculating the source position 
for each point on the grid. Every set of three points in 
the image plane defines a triangle and corresponds to a 
triangle in the source plane. The ratio of the area of the 
image plane triangle to the source plane triangle defines 
the magnification of that area of the source plane. 

In Figs.[l]&|2] we show the magnification maps for four 
realizations of minimum and four realizations of saddle 
points. Each map shows a small patch of the source 
plane, the color corresponding to the total magnifica- 
tion of the image (including any micro- or nanoimages) 
for a source at that position. Given the small size of 
the patches and the stochasticity of the placement of the 



4 



Chen & Koushiappas 



stars and halos, different realizations with the same prop- 
erties look significantly different. 

As previous works have also shown (e.g., 'Schechter 
fc Wambsganss|[2QQ2[ ), the magnification maps show a 
wealth ot features that can be described as of two types: 
peaks (or cusps) of high magnification, and valleys of 
lower magnification. In particular, the saddle points - 
which are preferentially demagnified by perturbations - 
show large valleys of very low magnification. In compari- 
son to the maps that do not contain subhalos, maps with 
subsolar mass halos show a significant excess of small 
structures, and the saddle point realizations have several 
smaller, even lower magnification, valleys layered on top 
of the larger valleys. 

4. MOCK LIGHT CURVES 

The magnification maps offer a detailed look into the 
lensing effects of subsolar mass dark matter halos. Ob- 
servations, however, are limited to sampling a light curve 
as the source moves across the magnification map. A sin- 
gle light curve will contain a small fraction of the infor- 
mation in the magnification map. Any distinct feature, 
such as a valley, that is observed in a track crossing in 
one direction will look very different if the source were to 
move in a different direction. In addition, finite sources 
will smear out some features. Finally, the observable is 
fluxes - and not magnifications - of images. Connecting 
observed fluxes to predicted magnifications require a de- 
tailed understanding of the smooth lens component, the 
lens environment, and the time delay between images. 

We create a catalog of mock light curves from the re- 
alizations of magnification maps. The effective source- 
plane velocity is 



A.. 



vi 



Do 



1 + zi Doi 1 



(3) 



where Vq^ vi^ and Vs are the velocities of the observer, 
the lens, and the source; zi and Zs are the redshifts of 
the lens and the source; and Doi^ Dis^ and Dqs are the 
angular diameter distances of the observer to the source, 
the lens to the source, and the observer to the source. 
We set a fiducial effective velocity of 600 km s~^. 

The source is modeled as a uniform disk, with fiducial 
rad ius of 2 x lO"*^^ cm. While the source profile is simplis- 
tic, [^b^ns^^FaL (2005) show that the microlensing 
signal is mostly sensitive to the half-light radius of the 
disk and the details of the disk profile are less important. 
The size of the quasar source depends on the wavelength 
of the observation. The smallest sources are observations 
of optical c ontinuum flux, ~ 10^^ cm (see, e.g., Morgan 
et al.||2010 '). Our fiducial disk size is smaller than previ- 
ously measured values but falls near the range of values 
expected and is larger than the radius of the last stable 
orbit for a Schwarzschild black hole of 2 x IO^Mq. 

Each magnification map is 1.5 x 10~^ arcsec on a side. 
At the fiducial velocity, an interval of 10 years will tra- 
verse an angular distance of 7 x 10~^ arcsec. Therefore, 
based on the size of the magnification maps, we can cre- 
ate light curves that contains 10 years of data sampled 
at week intervals. During that time interval the source 
moves in a trajectory that is randomly oriented with re- 
spect to the magnification map. 

We create 200 such light curves for each realization. 
Figure [3] shows sample light curves. Here some typical 



features are readily apparent: long planes of low magni- 
fication and cusps of high magnification at caustic cross- 
ings. In addition, some features can be identified as 
caused by subsolar mass halos: a small valley at week 
~ 350 in the upper panel and several small peaks be- 
tween week 175 and week 375 in the lower panel. 

Histograms of the magnification maps are shown in 
Figure |4] Any differences in the histograms between re- 
alizations with subsolar mass halos from those without is 
small: a slight excess of lower magnifications both saddle 
and minimum points, a slight excess of higher magnifica- 
tions in the minimum points. Without knowing the mag- 
nifications (and not the observed fluxes) of the images, 
these small differences may be impossible to measure at 
all. It is clear, then, that a simple statistical analysis of 
the measured fluxes in observed light curves will be sim- 
ilarly insufficient in detecting subsolar mass DM halos. 
We must, instead, seek an analysis of the features of the 
light curves in order to find more significant differences 
between realizations with dark matter halos and those 
without. 

5. ANALYZING LIGHT CURVE FEATURES 

Our mock light curves are without observational errors, 
and so a simple analysis can be utilized. We find the 
centers of peaks and valleys by the derivative of the light 
curve. The beginning and the end of a peak or valley 
is found by the curve's inflection points. Events that 
extend beyond the length of the light curve are ignored. 
Peaks and valleys which are not at least A mag = 0.01 
are discarded. A mag = 0.01 is significantly larger than 
the precision of the lensing simulation and similar to that 
expected from future observations. Accounting for all 
peak and valley events, some fraction of each light curve 
will lie between events. These areas will be referred to as 
'plains,' and they may be asymptotically rising or falling 
or flat. Figure [3] labels the peaks and valleys on a sample 
light curve. 

Table [2] shows the results of the analysis. The num- 
ber of events is similar for peaks, valleys, and plains, 
although the realizations with dark matter halos have 
~50-100% more events than those without. 

A more in depth look into the characteristics of 
the events considers the magnitude displacement and 
timescale of events. A peak or valley event consists of 
a rising side, an extremum, and a falling side. We mea- 
sure two magnitude amplitudes: one between the start 
of an event and the extremum and one between end of 
the event and the extremum. We sort the magnitude 
amplitudes by size - each peak has a minimum mag- 
nitude difference and a maximum magnitude difference. 
For plains, the magnitude amplitude is measured as the 
difference between the start and the end of the event. 
The timescales are measured as the total duration of an 
event. 

Results for the magnitude amplitude and timescale of 
events is shown in Figures [5) [6) and (T] The peaks show 
that the magnitude amplitude of the event is smaller for 
the realizations with dark matter halos than those with- 
out, although the timescales are similar. Similar results 
are found in the plains. A much larger difference, how- 
ever, is seen in the valleys. Here the minimum magni- 
tude difference is much smaller for the realizations with 
dark matter halos than for those without. A milder ef- 
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x(10' arcsec) 



x(10' arcsec) 



Realizations 


Peaks 


Valleys 


Plains 


S1-S7 


5621 


4332 


5774 


Shl-Sh7 


7774 


6673 


8445 


M1-M5 


3170 


2130 


3639 


Mhl-Mh5 


6027 


5027 


6268 



TABLE 2 
Event Statistics 



Fig. 1. — Sample magnification maps for a saddle point image. The realizations on the left (Shi, top; Sh2, bottom) contain subsolar mass 
halos, while those on the right (SI, top; S2, bottom) are without. Note that the realizations to the left have more fine structure than the 
realizations to the right. 

els of Figure [Sj Here the minimum magnitude amph- 
tude and the timescale for vaheys are plotted in contours. 
The realizations with and without dark matter halos are 
found in overlapping regions. This is unsurprising, since 
a large fraction of the the events in any realization is 
due to stars. As we consider smaller fractions of the 
events, we see that 50% of the events are found in dif- 
ferent parts of the parameter space. Events due to dark 
matter halos are found with minimum magnitude ampli- 
tudes less than 0.1 mag and event timescales less than a 
few months. Events due solely to stars display a range 
of minimum magnitude differences and timescales. This 
suggests that it would be impossible to classify a single 
event as being produced by dark matter halos or by stars. 
At the 50% level, though, valleys in realizations without 
dark matter halos have a minimum magnitude difference 
~ 0.3 and timescales on order of years. 



feet is seen in the maximum magnitude difference and in 
the timescale; the magnitude difference and the timescale 
are somewhat smaller in the case with dark matter halos 
than in the case without. Here it seems that when dis- 
tinguishing between microlensing events and events due 
to subsolar mass dark matter halos, valley events will be 
more useful. 

These differences can be summarized in the top pan- 
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x(10' arcsec) 



X (10 arcsec) 



Fig. 2. — Sample magnificatixon maps for a minimum point image. The realizations on the left (Mhl, top; Mh2, bottom) contain subsolar 
mass halos, while those on the right (Ml, top; M2, bottom) are without. Note that the realizations to the left have more fine structure than 
the realizations to the right. 



Figure |8] also shows the results of varying some of the 
fiducial conditions, such as source size, observation fre- 
quency, and source velocity. As the source size increases, 
events are smoothed out. Fewer total events will be 
measured, and the remaining events have smaller mag- 
nitude amplitudes. A distinct signal from dark matter 
halo events will be more difficult to discern. A less fre- 
quent monitoring of a lensing system will also obscure the 
differences between stars and dark matter halos; obser- 
vations will be sensitive to only longer duration events. 
A faster effective source velocity (900 km/s) may have 
similar effects to less frequent observations - although, 
in this case, the effect is mild - but it will sample from 
a larger area in the source plane (and more events) than 
the fiducial velocity can. 

6. DISCUSSION 



The light curve analysis shows that subsolar mass dark 
matter halos in galaxies leave clearly observable effects. 
They create 'nanolensing' events that are seen as small 
magnitude amplitude and short duration events in a light 
curve. The amplitude in magnitude is a less than A mag 
= 0.1 and require observations with small photometric 
errors. The timescales are significantly less than 1 year 
and require a short cadence of observations. 

There remains some challenges that may impair our 
ability to detect subsolar mass dark matter halos with 
gravitational lensing such as the dark matter halo profile, 
the fraction of surface density in stars and in subsolar 
mass halos, and the intrinsic variability of quasar sources. 

As discussed previously, our fiducial model uses H oc 
R~^-^ which is steeper than the p ex r~^"^ that most re- 
cent simulations suggest, although the possibility still re- 
mains that these dark matter halos will have higher cen- 
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Fig. 3. — Sample 10-year light curve with observations taken 
weekly drawn from Sh2 (top) and from Mhl (bottom). Peaks are 
labeled in red points and valleys are labeled in blue points. The 
extrema for those events are labeled in cyan and magenta points, 
respectively. Portions of the light curves remain unlabeled; they 
are referred to as 'plains,' although they may be asysmptotically 
rising or falling. The effects of subolar mass halos can be seen as 
a small valley around week ~ 350 in the top panel and the three 
small peaks between weeks 175 and 375 in the bottom panel. 
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Fig. 4. — Magnification histogram for all saddle point realiza- 
tions (left) and all minimum point realizations (right). The filled 
red bars represent the realizations with subsolar mass halos, while 
the black lines represent the realizations without. Simple magni- 
fication histograms are unable to distinguish between realizations 
with subsolar mass halos and those without. 



left uses the fiducial S ex R~^-^^ while the one on the 
right uses S oc R~^. These maps are lower resolution 
and show a larger area than in the realizations used in 
the light curve analysis, but it is clear that a shallower 
density profile will result in smaller structures, which will 
be more difficult to detect. 

While the density profile of subsolar mass halos has a 
significant effect on the size of features at ffxed mass func- 
tion and surface mass density, the same kinds of features 
in magniff cation maps - peaks, valleys, plains - appear 
regardless of the profile. In addition, any given magnifi- 
cation map is not a unique representation of the profile 
of the subsolar mass halos. For example, the two panels 
in Figure |9] could be made similar by increasing the to- 
tal mass or the subsolar mass halos in right panel. This 
degeneracy is not confined to extended halos either. The 
magnification map created by a distribution of fiducial 
subsolar mass halos can be approximated by a distribu- 
tion of subsolar mass point-like objects with a smaller 
surface mass density. 

The fraction of surface density in stars and in dark 
matter halos will also impact the observability of subso- 
lar mass dark matter halos. As the fraction of surface 
density in stars increases, there will be a smaller amount 
of mass that can be attributed to subsolar mass dark 
matter halos. If the stellar fraction decreases, more mass 
can be attributed to dark matter halos. The microlens- 
ing ffuctuations, however, will remain signiffcant until 
the stellar fraction drops below ~10%. 

Ignoring the effect of stars, as we vary the surface mass 
density in dark matter between a smooth component and 
a component of dark matter halos, we will see simila r 
effects as described by Schechter & Wambsganss ( 2002 ) , 
with only the size of the events on the source plane scaled 
down. As the fraction in dark matter halos grows, a 
growing number peak events interlace, crowding out the 
deep valleys seen in the saddle point magnification maps. 
Given that subsolar mass halos are not point sources, this 
effect w ill never been as strong as in the case described by 
[Schechter & Wambsganss (2002). As the fraction of dark 
matter halos decreases toward 0%, the peaks become less 
frequent, more distinct, and separated by large plains or 
valleys. 

An additional source of variability in observed light 
curves is the intrinsic variation of the quasar source. In 
observed lenses, light curves from different images of the 
same source are generally differenced in order to elimi- 
nate the effect of intrinsic variability. However, this ob- 
scures the differences between peak and valley events due 
to lensing. We can compare the variability due to lens- 
ing to the intrinsic variability by studying a sample of 
unlensed quasars. Schmi dt et al., (2010 ) characterize a 
sample of quasars taken from the Sloan Digital Sky Sur- 
vey (SDSS) Stripe 82. For each object, they quantify the 
mean variability amplitude as a function of the time be- 
tween observati ons by calculating th e structure function 
dCristiani et al.| |1996| [Giveon et al]|1999D , 



tral densities than calculated ( Ishiyama et al.|20"TQ ). The 
strong 'nanolensing' effects or subsolar mass dark mat- 
ter halos are dependent upon the density profile. Figure 
|9] shows magnification maps in which the entire surface 
mass density is put into subsolar mass dark matter ha- 
los using the fiducial mass function. The panel on the 



V{At) 



(4) 



At 



where the average Qai is taken over all pairs of observa- 
tions i and .7, whose time difference falls into the bin At. 



Schmidt et al.| ( fiOlO) model the shape of the structure 
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Fig. 5. — Statistics for 'peaks' in saddle point realizations (top) and minimum point realizations (bottom). Each event consists of a rising 
side, a maximum, and a falling side. Left: The smaller change in magnitude over the event (either the rising side or the falling side). 
Center: The larger change in magnitude. Right: The duration of the event in weeks. The filled red bars represent the realizations with 
subsolar mass halos, while the black lines represent the realizations without. 



function for a sample of quasars as a power law with am- 
plitude A and slope 7, Vmodei — A(At/lyear)^ and find 
that 0.07 < A < 0.25 and 0.15 < 7 < 0.5 contain over 
90% of the quasars. 

We calculate the structure function for light curves 
with and without subsolar mass dark matter halos. The 
results for the mean and the standard deviation are 



shown in Figure 10 with comparison lines for the power 
laws with the largest amplitudes, A = 0.25 and the 
bounds of the slopes 7 = 0.15 and 7 = 0.5. It is clear 
that the structure function is not particularly sensitive 
to the differences between lensing due to stars and lens- 
ing due to subsolar mass dark matter halos. We can 
see, though, that lensing causes a larger amplitude in 
variability compared to intrinsic variability at scales of 
> 0.3 years. At the smallest scales of < 0.1 years, on the 
other hand, the signal from intrinsic variability will dwarf 
that of lensing and image differencing will be necessary. 
At at intermediate scales - where we may hope to find 
nanolensing events - it is unclear if intrinsic variability 
could be confused with lensing events. 

We create a sample of color-selected quasars and a 
comparison sample of c olor-selected F/G stars using cri- 
teria based on those in "Schmidt et al." ('2010') using the 
SPSS ([Ab aza.iian et al. 2009; York et al. 2000; Gunn 
et al.|2005| |1998f ^Fukugita et al. 1996; Hogg et al. 2001). 
I'he cadence of observations of Stripe 82 is such that each 
object is sampled for 2-3 months per year (the median 
time is ~10 weeks). This sampling may be long enough 
to search the quasar light curves for peaks and valleys in 
a similar manner to that employed on the mock lensed 
light curves. Unlike mock observations, however, the r- 
band observations have typical magnitude errors of 0.03 
and the ^-band of 0.04. We create a simple peak and val- 
ley finder that looks for extrema in the r-band with dura- 



tions longer than 1 observation and amplitudes above the 
observed photometric errors and that have a correspond- 
ing event in the ^-band data. The quasar sample finds 
events with median duration of '^5 weeks, while the star 
sample has a median duration of ~3 weeks. The duration 
of quasar events is typically less than the total amount 
of time sampled, less than the typical duration of micro- 
and nanolensing events, and greater than what we may 
presume is noise found in the non- varying star sample. 
This could be an indication that events created by intrin- 
sic variability have different properties than those caused 
by lensing. Without a more richly sampled dataset with 
smaller photometric errors, however, any definitive con- 
clusions about the nature of intrinsic quasar variability 
is not possible. 

7. CONCLUSIONS 

Subsolar mass dark matter halos are a prediction of 
cold dark matter cosmologies and a test thereof. In ad- 
dition, measuring the small scale structure of dark mat- 
ter will aid in understanding structure formation and in 
making predictions for current and future dark matter 
particle detection experiments. 

Here we suggest that strong gravitational lensing could 
be used to detect subsolar mass dark matter halos and 
show that such halos create observable lensing effects. 
We describe how the effect of subsolar mass dark matter 
halos differ from the effect of stars in the galaxy lens: 

1. Subsolar mass dark matter halos add lensing events 
to light curves that are already populated by events 
from stars. 

2. Dark matter halos cause events with preferentially 
smaller amplitudes and timescales, particularly in 
the case of 'valley,' or dimming, events. In such 
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Fig. 6. — Statistics for 'valleys' in saddle point realizations (top) and minimum point realizations (bottom). Each event consists of a 
falling side, a minimum, and a rising side. Left: The smaller change in magnitude over the event (either the rising side or the falling side). 
Center: The larger change in magnitude. Right: The duration of the event in weeks. The filled red bars represent the realizations with 
subsolar mass halos, while the black lines represent the realizations without. 

curves: light curves from two images are offset by the 
measured time delay and differenced to eliminate the ef- 
fect of intrinsic source variability. Thus, a brightening 
event attributed to one image of a pair can be inter- 
preted as a dimming event in the other image. And a 
single event of ~ 100 days could be a peak event, a valley 
^ event, or some linear combination of both, and it could 
be due to stars or due to subsolar mass dark matter ha- 
los. 

Monit oring of Q2237+ 305 s howed a single >100 day 




10 100 

Time (weeks) 



10 

Time (weeks) 



Fig. 7. — Statistics for 'plains' in saddle point realizations (top) 
and minimum point realizations (bottom). Left: The change in 
magnitude over the event. Right: The duration of the event in 
weeks. The filled red bars represent the realizations with subsolar 
mass halos, while the black lines represent the realizations without. 

events, the minimum magnitude amplitudes are 
typically less than 0.1 mag and the timescales less 
than ^100 days. 

These effects are visible with a simple analysis of light 
curves and without knowledge about the larger-scale con- 
figuration of the lens that would require detailed lens 
modeling, such as the shape and size of the galaxy lens 
in which the subhalos are embedded or the intrinsic flux 
of the source. 

Several observed lens systems have shown variability 
at timescales less than '^100 days. These events are of- 
ten identified as 'brightening' events in differential light 



event in [Schmidt et al.| ([2002|. Q0957-h561 has been 
studied by several groups and at different times who each 
found a single rapid event wit h durations ranging from 
less than 1 day to tens of da ys ( iSchmidt fc Wambsganssl 
[1998; Gil-Merino et al. 2001l (Colley & Schild 2003). SBS 
1^20+530 and HE 1104-1805 have both been observed to 
have multip le events with durations of ^5 day s ( ^Burud) 
et al.||2002^ and -1 month ( |Ofek fc Maoz||2003D , respec- 



tively. it could be that none of these events are due to 
lensing by subsolar mass dark matter halos, and larger 
samples of monitored lens systems will be necessary in 
order to detect the effects of subsolar mass dark matter 
halos with any certainty. 

Future surveys may detect large number of quasars 
lensed by galaxies with sufficient monitoring to mea- 
sure precise time delays. The Panoramic Survey Tele- 
scope & Rapid Response System (Pan-STARRSf] and 
the Large Synoptic Survey Telescope (LSStJ^ are ex- 
pected to detect ^2000 and ~8000 lenses, respectively 
( |OgurifcMarshall|2010D . Both Pan-STARRS and LSST 
will perform monitoring and may measure time delays 
with an accuracy of a few days. A proposed dedicated 



^ http:/ /pan-starrs. ifa.hawaii.edu. 
^ http://www.lsst.org 
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Fig. 8. — The minimum magnitude amplitude - timescale plane 
for minimum source points (left) and for saddle point (right). The 
realizations with dark matter halos are shown in greyscale and the 
realizations without are shown in contour lines. In both cases, 
the contours show where 90%, 70% , 50%, and 30% of the events 
lie. From the top, the rows show the fiducial case, the case for 
which the source radius is 8 times larger, the case for which the 
observation frequency is 6 weeks, and the case for which the source 
velocity is 900 km/s. 



lens monitoring mission, Observatory for Multi-Epoch 
Gravitational Lens Astrophysics (OMEGA; [Moustakas 



of ~0.1 days for 100 lenses. It is important to note that 
in conjunction with these datasets, new methods of anal- 
ysis need to be developed, implementing subsolar mass 
dark matter halos into Bayesian methods of constraining 
the characteristics of lens systems. 

In summary, we present a methodology that can be 
used to extract the gravitational nanolensing signal that 
may originate from the presence of subsolar mass dark 
matter halos in lensing galaxies. We show that such sig- 
nal may potentially be observable, and that the nanolens- 
ing effects can be distinguishable from the microlensing 
effects caused by stars. Future surveys will observe a 
sufficient number of quasar lenses that could be used to 
undertake a systematic search for the nanolensing sig- 
nal due to subsolar mass halos. A potential detection of 
these effects will be invaluable as it will shed light in the 
nature of dark matter, it will provide information about 
very early structure formation and will expand our un- 
derstanding of the energetic sources that power active 
galactic nuclei. 
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